ON GOWDY VACUUM SPACETIMES 

HANS RINGSTROM 



Abstract. By Fuchsian techniques, a large family of Gowdy vacuum space- 
times have been constructed for which one has detailed control over the as- 
ymptotic behaviour. In this paper we formulate a condition on initial data 
yielding the same form of asymptotics. 



1. Introduction 

This paper is concerned with the study of cosmological singularities. By a cosmo- 
logical spacetime we mean a globally hyperbolic Lorentz manifold with compact 
spatial Cauchy surfaces satisfying Einstein's equations. A singularity is charac- 
terized by causal geodesic incompleteness (assuming the spacetime satisfies some 
natural maximality condition). Causal geodesic incompleteness, and thus singu- 
larities, is guaranteed in general situations by the singularity theorems. However, 
the question of curvature blow up at the singularity, and the related question of 
strong cosmic censorship are a separate issue. The desire to answer these questions 
motivated this paper. 

Most of the work in the area of cosmological singularities has concerned the spatially 
homogeneous case. However, some classes of spatially inhomogeneous spacetimes 
have been studied analytically and numerically. In particular, the so called Gowdy 
spacetimes have received considerable attention. The reason for this is probably 
the fact that analyzing the Gowdy spacetimes is on the borderline of what is doable 
and what is not. These spacetimes were first introduced in [|| (see also Q), and 
in the basic questions concerning global existence were answered. We will take 



the Gowdy vacuum spacetimes on R x T 3 to be metrics of the form (1.1). However, 
some sort of motivation for this choice seems to be in order. Below, we give a rough 
description of more natural conditions that lead to this form of metric. In fact, 



the conditions below do not imply the form (1.1), see || pp. 116-117. However 



the discrepancy can be eliminated by a coordinate transformation which is local in 
space. Combining this observation with domain of dependence arguments hopefully 
conv inces the reader that nothing essential is lost by considering metrics of the form 
( fi~l| ) . The description below is brief and we refer the interested reader to || and [|j 
for more details. The following conditions can be used to define the Gowdy vacuum 
spacetimes: 



It is an orientable globally hyperbolic vacuum spacetime. 



• It has compact spatial Cauchy surfaces. 

• There is a smooth effective group action of U(l) x £7(1) on the Cauchy surfaces 
under which the metric is invariant. 

• The twist constants vanish. 



2 



HANS RINGSTROM 



Let us explain the terminology. A group action of a Lie group G on a manifold 
M is effective if gp — p for all p £ M implies g = e. Due to the existence of 
the symmetries we get two Killing fields. Let us call them X and Y. The twist 
constants arc defined by 

k x = e a p 7S X a Y p VX s and k y = e^sX^^'Y 5 . 

The fact that they are constants is due to the field equations. By the existence 
of the effective group action, one can draw the conclusion that the spatial Cauchy 
surfaces have topology T 3 , S 3 , S 2 x S 1 or a Lens space. In all the cases except 
T 3 , the twist constants have to vanish. However, in the case of T 3 this need not 
be true, and the condition that they vanish is the most unnatural of the ones on 
the list above. There is however a reason for separating the two cases. Considering 
the case of T 3 spatial Cauchy surfaces, numerical studies indicate that the Gowdy 
case is convergent Q and the general case is oscillatory Q . Analytically analyzing 
the case with non-zero twist constants can therefore reasonably be expected to be 
significantly more difficult than the Gowdy case. In this paper we will consider the 
Gowdy T 3 case. 

Due to the numerical studies, cf. 0, the picture as to what should happen is quite 
clear. In order to formulate the conclusions, we need to parametrize the metric. 
One way of doing so is 

(1.1) 

g = e {T - x)/2 {-e- 2T dT 2 + d6 2 ) + e- T [e p da 2 + 2e p Qdadd + {e p Q 2 + e - p )dS 2 }. 
Here, r £ K and (6, a, 6) are coordinates on T 3 . The evolution equations become 

(1.2) Prr-e- 2T P ee -e 2P (Q 2 T ~e- 2T Q 2 ) = 

(1.3) Q TT -e- 2T Q 9e + 2(P T Q T ~e- 2T PeQe) = 0, 
and the constraints 

(1.4) A T = P 2 + e- 2T P 2 + e 2P (Q 2 T + e- 2T Q 2 ) 

(1.5) \ B = 2(P P T + e 2P Q e Q T ). 

Obviously, the constraints are decoupled from the evolution equations, excepting 



the condition on P and Q implied by (1.5). T hus the equations of interest are 



the two non- linear coupled wave equations (L2)-(1.3). In this parametrization, the 



singularity corresponds to r — > oo, and the subject of this article is the asymptotics 



of solutions to (1.2)- (1.3) as r — > oo. The asymptotics we derive will then be used 
to obtain conclusions concerning curvature blow up. There is a special case of 
these equations determined by the condition Q = which is called the polarized 
case. This has been handled in |j, which also considers the other topologies. The 
asymptotic behaviour of the solution P to (fT^) in the polarized case is given by 

P(r, 6) = v(6)t + (f>{9) + e~ eT u(T, 6) 

where e > and u(r, 6) — > as r — > 0. In this situation, v and 4> are arbitrary 
smooth functions on the circle. In the general case, the numerical studies indicate 
that the "velocity" v should typically be confined to the open interval (0, 1). To be 
more precise, the following asymptotics are expected in general: 

(1.6) P(t,6) = v(6)T + (f>(9)+e- eT u(6,T) 

(1.7) Q(t,9) = q(0)+e- 2v ^m9)+Mr,O)} 
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where e > and w, u — > as r — > oo and < v(6) < 1. A heuristic argument 
motivating the condition on the velocity can be found in . However, the numerical 
simulations also indicate the occurrences of "spikes". Let us describe one sort of 
spike that can occur. It can happen that at a spatial point 8q, P t will, in the 
limit, have a value greater than 1 whereas the limiting values of P T in a punctured 
neighbourhood of 8q, P t belong to (0, 1). Furthermore, Q converges nicely in a 
neighbourhood of 8q with a zero of the spatial derivative of Q at 6*o ■ Beyond the 
numerical indications of these types of features, a family of solutions with spikes 
have been constructed in Jj"l| ], so that the behaviour described above is known to 
occur. The type of spike described above is a "true" spike. There are also other 
types of spikes called "false" spikes at which Q has a discontinuity. We refer the 
reader to |llj for more details. One relevant question to ask is whether the spikes are 
a result of a bad parametrization of the metric or if they really have a geometrical 
significance. It seems that the false spikes are a result of bad parametrization, but 
the true spikes can be detected by curvature invariants 

In this paper we will not be concerned with spikes, but will focus on solutions 
with an asymptotic behaviour of the form (1.6)- (1.7). By the so called Fuchsian 
techniques one can construct a large family of solutions with such asymptotic be- 
haviour. In fact, given functions v, 4>, q and ip from S 1 to R of a suitable degree of 
s mo ot hne ss and subject to the condition < v < 1 , one can construct solutions to 
(F2)-(L3) with asymptotics of the form (1.6)-(1.7). The proof of this in the real 
analytic case can be found in |^| and |[o| covers the smooth case. One nice feature 
of this construction is the fact that one gets to s peci fy fou r functions freely, just as 
as if though one were specifying initial data for (1.2)-(1.3). 

The purpose of this pa per is to p rovide a condition on the initial data yielding the 
asymptotic behaviour (l.€)-(1.7). There are several reasons for wanting to prove 
such a statement. As was mentioned above, one can construct a large family of 
solutions with the desired asymptotic behaviour, but it is not clear how big this 
family is in terms of initial data. In this paper we prove the existence of an open 
set of init ial data yielding the desired asymptotics. Observe that the equations 
( |l.2| )-( |L~3| ) are not time translation invariant, so that at which time one starts is of 
relevance. The open set in the initial data will thus depend on the starting time 
tq. The condition demanded in this paper is not only sufficient, but in fact also 
necessary for obtaining asymptotics of the form ( 1.6)-( 1.7) in the sense that, if a 
solution has this form of asymptotics, then for a late enough time, the condition on 
the initial data will be satisfied. In this sense, the condition desc ribe d in this article 
is a characterization of the solutions with asymptotic behaviour ([L6])-([L7]) in terms 
of initial data. Observe finally that the condition, even though it is formulated as a 
global condition on all of S 1 in this paper, can be applied locally due to domain of 
dependence arguments. Thus the condition prescribed here should be of relevance, 
and should in fact be applicable in a neighbourhood of almost all spatial points, 
even in the case with spikes. The problem in the general case of course being that 
of proving that the evolution takes you to such a region. Finally, let us observe that 
this is not the first result in this direction. In |4j , Chrusciel considers developments 
of perturbations of initial data for the Kasner (|, |, — |) metric within the Gowdy 
class. He proves, among other things, curvature blow up for these developments. 
In our setting the Kasner (|, |, — i) metrics correspond to P = Q = 0. 
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2. The equations and an outline of the argument 

Since the method does not depend on the dimension, and since the arguments in 
a sense become more transparent in a more general setting, we will consider the 
following equations on R x T d : 

(2.1) P TT -e- 2T AP-e 2P (Q2 _ e - 2r |VQ| 2 ) = 

P{to,-) = Pa, Pt(t ,-) = Pi 

(2.2) Q TT - e- 2T AQ + 2{P T Q T - e~ 2T VP ■ VQ) = 

Q{to,') = 9o, Qt(to,-) = Qi- 
These equations have some similarities with wave maps. Let 

d 

g = -dt (g> dt + ^^dx l ® dx % 
i=i 

be the Minkowski metric on R x T d and let 

go = dP®dP + e 2P dQ dQ 

be a metric on M 2 . Observe that (K 2 , go) is isometric to hyperbolic space. The wave 
map equations for a map from (fx T , g) to (R 2 , go) is given by the Euler-Lagrange 
equations corresponding to the action 

J go,a b g a(3 d a u a dpu b dtdx = J[-P 2 + |VP| 2 ~ e 2P Q 2 + e 2P \VQ\ 2 ]dxdt. 

This should be compared with (2.1) and (2.2) which are obtained as the Euler- 
Lagrange equations corresponding to the action 

J [-P 2 + e- 2r |VP| 2 - e 2P Q 2 T + e 2p - 2T \VQ\ 2 }dxdr. 

The exact statement of the result can be found in section ^. It is however rather 
lengthy, and therefore we wish to state a somewhat less technical consequence 
here. We need to define some energy norms, but first note the following convention 
concerning multi-indices. If a = (ai, ay) where on are non-negative integers, 
then 

q\<A f 

D a f = — t 

J dxl 1 ■ ■ ■ dx a d d ' 

where \a\ = oti + ... + ad- The natural energy norm for P is 

(2.3) £ k (r) =£ k (P :T ) = lj2 I l(D a d T P) 2 + | V D a P\ 2 ]d6 

1 \ a \=k jT " 

and the one for Q is 

(2.4) E k {r)=E k {P 1 Q, T )= l - V / [e 2P (D a d T Q) 2 + e 2p - 2T \VD a Q\ 2 ]de. 

\Oi\ — k 

Finally it is natural to introduce 



T k {r)=T k (P,r)= l - I |V^ Q P| 



dQ. 
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If Po,Pi,qo,qi £ C°° (T d , R) , we will by e k (po ,Pi,t) mean £ k with P replaced by po 
and P T replaced by p\. We associate e k {po, qo, Qx, t) with E k and v k (po) with jFfc 
similarly. 

Theorem 2.1. Letp , Pl £ C°°(T d ,R) sate/y 

27 < Pi < 1 - 2 7 , 

where 7 > 0. Then, if To is big enough and e k (po, go , qi, tq), fc — 0, = 2[d/2] + 



2 are small enough, we have smooth solutions to l[2.\ ) and (2.i) on [ro,oo) x T d 
with the following properties: there are v,w,q £ C°° (T d ,R) with 

7 < v < 1 — 7 

and polynomials iti.k, 7r 2 .fc in r — To /or every non-negative integer k such that 
\\P - p|lc fc (T d ,K) < Ti.fc exp[-27(r - r )], 

^ || e 2 "[-D a Q - D a q}\\ c(Tdm < n 2 , k , 

\a\<k 

where p = v ■ (r — t ) + w. 

Remark. The sizes of To and the e k (po, qo, Qi, To) on ly depend on 7, efc(pojf>ij To) 
and Vkipa), k — 0, ...,md- The reader interested in more conclusions is referred to 
section 0. 



Let us give an outline of the proof. Given a smooth function P, we can view (2.2) 
as a linear equation for Q. The central part of the argument is an analysis of this 
linear equation assuming P satisfies 

0<7<P r <l-7 

and £k(P) bounded for k — 0, ...,rrid — 2[d/2] + 2. Under these conditions on P 
one can prove that E k satisfies 

(2.5) El /2 {r) < Vu{t - T )exp[- 7 (T - t )] 

for k = 0, ...,to^. Here V k is a polynomial in t — To whose coefficients depend on 
the values of E kl £ k and T k at r = to for k = 0, ...,nid- The polynomials P^ have 
one important property. If one lets £7/. (to) go to zero, then the coefficients of the 
polynomial go to zero, assuming the other parameters are constant. In consequence, 
under these conditions on P, one has very good control of Q. The idea is then to 
consider the following iteration: 

(2.6) P n>TT -e- 2T AP n -e 2P —(Ql T -e- 2T \VQ n \ 2 ) = 

(2.7) Q„, rT -e- 2T AQ„ + 2(P„_ 1 , T Q n , T -e- 2r VP„_ 1 -VQ n ) = 

(2.8) P , rr - e- 2r AP = 
where 

P«(t , ') = PO, Qn(r0, ■) = qQ, Prur(T" , •) = Pi, Qn,r(T , •) = 

Observe that we only need to solve linear equations during the iteration in that 



if P„_i is given, then we can compute Q n using (2.7), so that (2.E) also becomes 
a linear equation. One sets up an induction hypothesis on P n amounting to the 
statements 

(2.9) <7 < P„ >T < I-7 and £ k {P n , t) < c k < 00, k = 0,...,m d . 
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It is not too difficult finding conditions on the initial data of P ensuring that Pq 
satisfies these conditions. By the above observations concerning the linear equation 
(p7^), one gets very good control of the behaviour of Q n+ i if ( fHi| ) holds. Insert- 
ing this information into ( |2.6| ), it turns out that the propagation of the inductive 
hypotheses essentially boils down to a smallness condition on E^ir®) k = 0, 



due to the structure (2.5). In this way, we produce a sequence of iterates obeying 
( [2.9| ). The arguments proving convergence turn out to be similar to the arguments 
proving the propagation of the inductive hypotheses, and a smallness condition on 



Ek(T~o) ensures the desired behaviour. Thus one produces a solution to (2.1) and 
( pT^ ) with certain properties. Due to the fact that one knows the solution to have 
these extra properties, one can show that it has the desired asymptotic behaviour. 

3. Local existence 

Let us here state the local existence result we will need. We are interested in 
equations of the form 

(3.1) DP(t,x) = F(t, x, P, P t , VP) 

P(*o,-) = Po, Pt(io,-)=Pi- 



Proposition 3.1. Consider the equation (3.1) where to G R_, (t,x) G R_ x T , 
(p ,Pi) G H k+1 (T d ,R l ) x H k (T d ,R l ), F is a smooth function and k > m d /2 = 
[d/2] + 1. Here, R_ = (-oo,0). Then there are T U T 2 G R_ with T x < t < T 2 
such that there is a unique solution solution of \3.\ ) in 

(3.2) C(I,H k+1 (T d ,R 1 )) n C 1 (/,iZ' fc (T <i ,R J )) 

where I = [Ti, T2]. Let T ma x be the supremum of the times T G M_ such that there 
is a unique solution to j[3.j ) in jS.^j for I — [to, T] and define T m i n similarly. Then 
either T max = 0, or 

SUP ||P||c 1 ([0,t]xT«',Mi) = °°- 

te[to,T max ) 

The statement for T m j n is similar. 

Remark. The fact that T d is compact makes some of the usual conditions on F 
unnecessary. 

The proof uses estimates of the form (6.4.5)' of Q adapted to the torus case. A 
similar result for k > m,d only requires Sobolev embedding and is sufficient for our 
purposes. We will later solve the non-linear problem by carrying out an iteration, 
and it will then be of interest to solve equations of the form 

Qu - AQ = GiQt + G 2 ■ VQ + G 3 

Q(*o,-) = 9o Qt(to,')=Qi 

where Gi, G 2l G3 G C°°(R_ x T d , R) and the initial data are smooth. Observe that 
local as well as global existence of smooth solutions to this problem is assured by 



Proposition |3.1| . Observe finally that equations of the form 

P TT (r,a;) - e - 2T AP(r, x) = F(t, x, P, P t , VP) 
on R x T d can be transformed to equations of the form (pTj) by the transformation 
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4. The polarized case 

It is always instructive to start with an easier subcase, and in this case we have the 
added incentive that the zeroth iterate of the iteration belongs to this subclass, so 
let us consider the polarized case. Let P solve 

(4.1) P TT - e- 2T AP = 0. 

Proposition 4.1. Consider a smooth solution P to the polarized equation ( [^. 
Then there are v,w £ C°°(T d ,M.) and positive constants C k £ K for every non- 
negative integer k such that 

(4.2) \\P T - w|| C fc( T d,R) + H-P -VT- w\\ C k (T d M) < C k (l + r)e~ 2T 
for all t > 0. 



Proof. The energy £ k defined by ( |2.3| ) satisfies 

(4.3) ^ < 0. 

ar 

By Sobolv embedding, we conclude that 

[\\D a d T P\\ c{T ^ m + e- T \\VD a P\\ c{TdM) ] < C k < oo 

\a\ = k 



for all k. Inserting this information into (4.1), we conclude that 

\\D a d 2 T P\\ c(Tdm < C k e- T 

\a\—k 

for all k and a. We conclude the existence of a function v G C°°(T d ,M.) such that 

HA- - v\\ C k( T d yM) < C k e' T , 

which in its turn proves the existence of a w £ C°°(T d ,R) such that 

\\P-vt- w\\ ck{TdM) < C k e~ T . 

One consequence of this is of course that 

\\P\\c*(T*M < C k {l + r) 

for r > 0. Inserting this in (4.1) and going through the same steps as above, one 
ends up with (fO|). □ 

Observe that in the end, it turns out that P and all its spatial derivatives do not 
grow faster than linearly. However, the natural consequence of the boundedness of 
£ k is that expressions of the form 



- V / e- 2T \VD a P\ 2 d9 
2 ^uJt-> 



2 

\a\ = k 

are bounded. In other words, the form of the energy ( |2.3[ ), forced upon us by the 
energy methods, is not well suited to the behaviour of the solutions. The way 
we achieved the linear growth estimate in the proposition above, was through a 
procedure which was very wasteful of derivatives. This is not likely to be successful 
in the general non-linear case. However, there is another point of view, and we will 
describe it in the next section. 
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5. Energies 

In this section we gather some observations concerning the type of energies we will 
be using. Let t G R + = [0,oo), I = [r ,oo) and let V G C°°(R x T d , K). Let 

(5.1) ^(V,r) = \Y. I VJD a 4>\ 2 d6. 

Lemma 5.1. Assume -0 G C°°(IR x T d ,R) satisfies 

£1 /2 (^,t) < e k < oo 



/or t £ I and fc = 0, m, where £ k is defined in (2.S ), m > rrid/2, md — 2[d/2] + 2 
and the e k are constants. Then, for k < m — rrid/2, 

[\\D a dMr, 0||o(r-« + e- T ||VZ>Xr, ■)||o(r*,R)] < C(e k + e k+md/2 ) 

\a\=k 

for t £ I, where the e k may be omitted if k > 0. Furthermore, 
(5.2) FI' 2 W,t) < C[^ fe 1/2 (^r )+e fe+1 (r-ro)] 

if k < m — 1, so that 

\\VD a 1p\\ C (T*,WL) < C [Fk+m d /2&> T 0) + £k+l+ md /2(T - T )] 

\a\—k 

if k < m — rrid/2 — 1. 

Remark. When we write ||/||c(T d .B) f° r a vector valued function /, we mean the 
sup norm of the Euclidean norm of the function. 

Proof. The first inequality follows from Sobolev embedding, as well as the third, 
given the second. The second is proved by computing 

= V / VD a i> ■ VD a d T il>d6 < 2CTl ,2 £ 1 k , + 1 . 

(XT J ' rpd 



\a\=k 



□ 



Observe that the main point of this lemma is the estimate fl5.2| ). We get a linear 
l li 

growth estimate for T k (ip, t) if we know that £k+i(ip, t) is bounded for r £ I. In 
other words, there is a price for this sort of estimate, but we only have to pay one 
derivative. 



As has already been mentioned, when considering (1.3), the following energy will 
be of interest, 

(5.3) %Cr)4E / [e 2 ^(D a d T 2 + e 2 ^\WD a e}dO. 

1 \a\=k JTd 

Lemma 5.2. Let r G R+ and I = [t 0j oo). Assume rj, ^ e C°°(l x T d , R) and 
that 

0<7<ry r <l-7<l 
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on I x T d where 7 is a constant. Then 

(5.4) ^M)<_2 7 £ fe +]T / UMD^dridO, 

T \ a \=k JTd 

where 

(5.5) f a ( V ,0 = d T (e 2r >D a d T - V • (e 2 "- 2T V£> Q 0- 

If a = a + ei, where e; is an element of Z d whose I :th component is 1 and whose 
remaining components are zero, we have the following recursion formula for f a : 

(5.6) f a = dif & - 2d lV f & - 2d l d T r 1 e 2 iD & d T £ + 2V(d lV ) ■ (e 2 ^ 2T WD & 0- 



Proof. Estimate 
dE k 
dr 



= Y] I [d T {e 2r) D a d T C)D a d T t-ri T e 2l i{D a d T C)' 2 



\a.\—k ' 

+ {n T - l)e 2 "- 2T |VL> Q £| 2 + e 2 "- 2T (VZ) Q ■ (V D a d T £,)]d6 < 
< -2jE k + V / [d T (e 2r <D a d T £) - V • (e 2r '- 2T VD a O}D a d T Zd8, 

it JT d 



and (5.4) follows. We have 

f Q = d T {e 2ri D a drC) - V • (e 2 "- 2r VL> a £) = 

= d T di{e 2n D & d T t) - d T (2d l r 1 e 2ri D & d T - «9/V ■ (e 2 ^- 2T VD & ^+ 

+V • {2d l ne 2 ^- 2T VD & ^) = dif & - 2d lV f & ~ 2d l d T ne 2 ^D & d T ^+ 

+2V{d l i 1 ) ■ (e 2, '- 2T VD & Z), 

proving (p^). 

6. Iteration 



Consider the iteration ( |2 . 6|) -(2.7 ) . We will only be interested in the future evolution 
of solutions to the corresponding non-linear partial differential equation, and we 
will implicitly assume the time interval on which our estimates are valid to be 
I = [to, 00). Let 

En P71 En — 1 j &n ,k — £k ( En : ' ) ; &n .k — £k ( E n , ■ ) , 

where £ k is defined by (2J3), and 

Qn = Qn~Qn-l, E n> k — E k (P n - 1 , Qn , ') , E n , k = E k (F n - 1 , Q n , •) 



where E k is defined in (5.3). Observe that these expressions are all independent of 
n if we evaluate them at To, that £ n .k is defined if n > 0, £ n ,k and E n .k arc defined 
if n > 1 and E n , k is defined if n > 2. 

Conditions and conventions concerning initial data. We only consider initial 
data (poiPij QotQi) with the property that there is a 7 > such that 

(6.1) < 2 7 < Pl (0) < 1 - 2 7 < 1 V 6 e T d . 
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Secondly, e k , v k and e k will be taken to be constants satisfying 
(6-2) S^iro) < e k , ^ k /2 (P n , r ) < u k , E^(r ) < e k 

for k = 0, nid, md — 2[d/2] + 2, with e k and e k positive. 
Induction hypothesis. We assume that 

(6.3) ^i 2 (r)< efc + lVre/ 
and that the following inequality is fulfilled on I x T d : 

(6.4) \P n>T {T,6) - P n , T {T ,Q)\ < A 7 . 

Here A7 should be suitably small relative to 7, but not too small. We will below 
assume A7 = 7/4 to hold. 

Observe that if ([O]) holds for n and m, then 

(6.5) |iV(r,0)-P m , T (T,0)| < 2A 7 . 



Observe also that (|6.4| ) implies that 

(6.6) < 2 7 - A7 < P n . T < 1- (2 7 - A7) < 1 

on I x T d . In the course of the argument, we will give inequalities involving 7, e k , v k 
and e k introduced in ( |6.l[ ) and ( |6.2| ) such that if they are fulfilled, the induction 
hypothesis is propagated. By imposing additional requirements, one can then prove 

ffi 71 f 1/2 < -f 1 / 2 

\ u -' ) n+l,max — °n,max 

where 

(6-8) = sup^fo + ■•• + SU P^m d - 

re; re/ 

Thus one obtains convergence. Observe that it makes sense to speak of the suprema 
once we have proven that ( |6.3| ) is valid for all n. Let us note some consequences of 
the induction hypothesis. 



Lemma 6.1. Assume (6.5) is satisfied up to and including n—l, k = 0, md and 
that (6.i) holds. Then 

(6.9) Yl l|V£ Q P„-i|U 2( T7 R ) <CK + (e fc +i + l)(r-ro)] 

\a.\ — k 

(6.10) J2 I|V£ q P„-i||l 2( t7m) < C£ l J\^ ■ (r - r ) 

|a|=fc 

on J, /or < k < md — 1 Finally, for n > 2, 

(6.11) |P B _i(r, 8)\ < C8Z\n^ ■ (t - to) 
on /. 

Proof . The estimates (6J5) and ( |3.1C| ) follow from Lemma 5d. In order to prove 
( |S.ll| ), we estimate 

\P n ^(r,9)-P n . 2 (r,9)\ = \ f P n __ ltT (a,0)da\ < C£ l J\ ma ^r - r ), 

where we have used Lemma 5.1 in order to obtain the last inequality. □ 
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Lemma 6.2. If (6.&) holds, then (6.Z) holds for n = 0. 



Proof See ((L3|). □ 

Lemma 6.3. There are constants Cd such that if (6A) and ((LS) are fulfilled for 
n — 1 and 



(6.12) 

if d = 1 , and 
(6.13) 

if d> 2, then 



on T a 



(e 2 + l)e- ro < ci 7 , 
[^ md /2+i + (e mi /2+2 + l)]e 

r OO 

/ e- 2T |AP„_i(r,-)|dr < A 7 /2 



~ 2ro < c dl 



Remark . Th e expression (1/2 + £3) exp(— 2ro) could also be used as the left hand 
side of ( 6.12 ) if one is prepared to keep track of one more derivative, c.f. the higher 
dimensional case of the argument presented in this paper. 

Proof. We get a division into two cases depending on the dimension d. If d > 2, 
then md > 4, so that md/2 + 1 < rrid — 1. In consequence, we can use the estimates 
of Lemma l^lj to obtain 

/•OO /"OO 

/ e~ 2T \AP n ^(T, -)\dr < C / e- 2 > md/2+1 + (e„ ld/2+2 + l)(r - r )]rfr < 

< CL/2+i + (e md /2+ 2 + l)]e- 2ro . 

If d = 1, we do not have enough control to ensure the linear growth of the third 
spatial derivative in the L 2 -norm, and therefore have to resort to using our control 
on £„_i,2- We get 

\\Pn-iftoWas 1 ) < C \\P n -ifiee\\ ^{s 1 ) < Ce T (e 2 + 1), 
where we have used ( ^3 ) . Thus 

/>oo 

/ e~ 2T \P n ^. ee (r,-)\dT < Ce- T "(e 2 + 1) 

J T 

and the lemma follows. □ 



Lemma 6.4. If the conditions of Lemma 6.c and (6.1 ) hold, the inductive hypothe- 
ses (6.c ) and \6.^J are satisfied for n = 0. 



Proof. The lemma follows by combining Lemma 3.2, 6.3 and ( 



7. The n:th step 



The first task is to estimate the behaviour of E n ,k- The point of the argument 
is to demand that all the iterates are such that P n , T belongs to a region where 
Q n ,T, at least intuitively, should decay to zero exponentially. Consequently, we 
hope to achieve an exponential decay for the energies E n ^- This is in fact the 
case, but as the argument is constructed, the natural estimate that appears is a 
polynomial times an exponentially decaying factor. The polynomials that appear 
have an important property we wish to formalize. 
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Definition 7.1. Let V be a polynomial in r — To depending on the e k and i/j 

for fc = 0, md. We say that V is a Q- dominated polynomial if the coefficients of 
V are polynomial in e k , and ^fc and go to zero when the go to zero while the 
other expressions are fixed. 

Remark. The Q-dominated polynomials we will consider will always be independent 
of n. 



Lemma 7.1. Assume that (6.&), (6.S) and (6.4) are satisfied for n — 1 and let 
< k < m c i. Then 

(7-1) E^(t) < V k (r - T ) exp[- 7 (T - t )] 

on I, where Vk is a Q-dominated polynomial. Furthermore, 

(7.2) Y, {\\^ Pn - 1Da 9rQn\\c(T^M) + ||e p »- 1 - T Vi? a Q n ||c(r^)} < 

\a\=k 

< Qfe+l(T-To)cxp[-7(T-T )] 

for k = 0, rrid/2, where the Qk are Q-dominated polynomials. 



Let us make some preliminary observations. By Lemma 5.2, we have 
dE n _k 



(7.3) 
where 



dr 



< -2 7 £ n , fc + V / fn, a D a d T Q n d0 



\a\—h 



fn,a = d T (e 2P '^D a d T Q n ) - V • (e 2P ^- 2T VD a Q n ). 

If a = a + ei, where e; is an element of Z d whose i:th component is 1 and whose 
remaining components are zero, we have the following recursion formula for f n . a - 

(7.4) / n , Q = dif nA - 2d l P n . 1 f nA - 2d l d T P n - 1 e 2P ^D & d r Q n + 

+2V(5 / P„_ 1 ) • (e 2P ^- 2T WD & Q n ). 

Lemma 7.2. If 

(7-5) ||e _P "- 1 /„,Q||L2 (T d !R) < TT k exp[-7(T - t )] 

for all a such that \a\ = k, where irk is a Q-dominated polynomial independent of 
n, then an inequality of the form (7.1) holds. 



Proof. By ( |73| ) and (7.5) we have 
dE n _k 



dr 



< -2 7 £„, fc + V / f n , a D a d T Q n d9 < 

I i i. *T d 



\ol\— k 



<-2-/E n , k + Y H e Pn ~ 1 fn,a\\L*(T*,n)\\e Pn ~ 1 D a d T Q n \\ L 2( T i i&) < 

\a\=k 



1/2 



Thus 



< -2 7 £ n , fe + V2Cir k exp[- 7 (T - t )]£, 



— {exp[2 7 (T - r )]E n<k } < V2C7r fc {cxp[2 7 (T - T )]£», fe } 1/2 
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which can be integrated to 

{exp[2 7 (r - T )]E ntk Y' 2 < ^ fc 2 (r ) + f 2- 1 / 2 Cir k (s - r Q )da 

and the lemma follows. 

Proof of Lemma \7. \ . Observe that 

fn,o = d T (e 2Pn - 1 d T Q r , 



„2P„ 



[Qn,rr ~t~ 2p 



- V • (e 2P - 



-It 



VQ r , 



2e~ 2T VP n 



VQ r , 



0. 



Let us also observe that when the two first terms in (7.4) hit an expression involving 
exp(2P„_i), then the effect is to differentiate the expression regarding the expo- 
nential expression mentioned as a constant. Inductively we thus get the conclusion 
that if \a\ = k + I, then /„ iQ consists of terms of the form 



(7.6) 

and 

(7.7) 



C Q;ij 



„2P„ 



■ 1 L> ai 9 r P I1 _ 1 L> Q2 9 T ( 



2P n . 



-It 



V£> ai P„_i ■ VD a H 



where |ai| > I and ai + a2 = a. We will have to use different estimates for different 
k:s. 



Zeroth order energy. Observe that 

oIt 



< -2 7 E n , 



so that 



1 /2 

E' (r) < e exp[- 7 (r - r ) 



for r > tq. Thus (7.1) holds for k — 0, with Vq = cq. 

Intermediate order energies. The condition I < k < m^/2 defines what we 
mean by intermediate energies. We carry out an argument by induction. Assume 
that we have (7.1) up to and including fc, < k < m^/2 — I. By Lemma 7.2, we 
need to consider 



,a\\L 2 {T d , 



for \a\ = k + 1. In order to deal with terms of the form (7.6), we need to estimate 



(7.8) 



e ^ D ^d T P n _ lD ^d T Q n \\ L ^ Td 



where a\-\-a2 = ol and |ai| > 1. By Lemma 5.1 and the induction hypothesis on n 
( |5.3| ) we can take out D ai d T P n -i in the sup norm. The remaining part is bounded 

— 1/2 

by y/2E n ', a ,. By the induction hypothesis on k and the fact that ja^l < k, we get 
the conclusion that 

\\ e P~-i D ^ dT P n _ lD «2 dT Q n \\ LHTdM) < 7rfc+iexp [„ 7 ( T _ To )]. 

The estimate for terms of the form ( [7.7p is similar, and we can thus apply Lemma 
7.2 in order to obtain (7.1) for the intermediate energies. 
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High order energies. By high order we mean md/2 + 1 < k < m^. Since we 
cannot assume to control the sup norm of the derivatives of P n -\ indefinitely, we 
need to change our method. Observe that if \a\ < k < md/2, then 



(7.9) 



\D a (e p '- 1 Q 



n,r )\\L 2 (T d TS 



< 7T 



)VQ n is 



for some polynomial 7r a , c.f. Lemma 5.1. The argument for exp(P„_i 
similar, and we get 

(7.10) ||e Pn - I <5„,x||c(T'i,R) + l|e P "" 1_T VQ I1 || c(Tci:K) < Q 1 exp[- 7 (r - r )] 

by S obolev embedding, where Q\ is a Q-dominated polynomial. Let us assume that 
( [7d| ) is satisfied up to and including rrid/2 + I < rrid — 1, and that (7.2) is satisfied 
for < k < I. For I — we know this to be true. We wish to prove an estim ate of 
the form (|7.5| ) and as for the intermediate energies, we need to consider ( |7.§| ) where 
Q'i + a2 = ot, \ot\ = + 1 + 1 and |ai| > 1. If |oti| < md/2, we can take out 

D ai d T P n -i in the sup norm to obtain 

\\ e P"-iD<*idrP n -lD aa drQn\\z*pW < ^2\\ d r P n . X \\ C (T*A) E n^ 2 | 

which satisfies a bound as in ([7^) by Lemma |5.l| , the induction hypotheses on I 
and the fact that \a%\ < md/2 + I. If \ct2\ < I we can take out e Pn - 1 D a2 d T Q n in 
the sup norm in order to achieve a similar bound using the induction hypothesis 



on I and the boundedness of 8. 



n,\ai\ ■ 



The argument for terms of the form (7.7) is 



similar. Since |ai| > md/2 and |aa| > I cannot occur at the same time, we are 
done. We have thus proven (7.1) for k = md/2 + 1 + 1. We now need to prove that 
fiT2\ ) holds for k = I + 1. However, this can be proven in the same way as (7.9); 
replace Q„ iT in that inequality with D a d T Q n for \a\ < I + 1. □ 

Let us now turn to the problem of estimating 



Lemma 7.3. Assume (6.i), (6.S) and (6.4) are fulfilled for n — 1 and let k — 
0, md- Then 

(7.11) £ l J 2 k {r) < £ l J 2 k {r ) + f V k (s - r ) exp[-2 7 ( S - r ))ds, 



where V k is a Q- dominated polynomial. 
Proof. Observe that 



dr 



It is thus of interest to estimate 



D<*{e^Qi ;T ) 



= E E C (3l ,...^ +2 D^P n _ 1 ---D^P n „ 1 e 2P ^D^d T Q n D^-d T Q n 

j=0 P 1 + ...+P j+2 =a 

in L 2 -norm. Consider a term in the sum. Observe that at most one \f3i | can be 
greater than md/2. Combining Lemma 6.1 and 7.1, we conclude that 

\\D a {e 2P ^Ql tT )\\ LHTdM) < 7r fc exp[-2 7 (r-ro)] 
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where n k is a Q-dominated polynomial. The argument for 

li 2 (T d ,R) 



\D a (e 2P ^- 2T \VQ n \ 2 ) 



is similar, and we obtain 



d£ n 



< 2V k exp[-2 7 (r - r )]£„ 



1/2 



fc ' 



which can be integrated to ( 7.11 



We now wish to specify conditions that imply ( |6.3| ) and ( p.4| ) for n. 

Lemma 7.4. There are constants Cfe, and non-negative integers iu and j k such that 
if 

(7.12) e k < c fc (l + eo + e md + v md )~ ik l jk 



for k = , TO,;, and i/ie relevant condition in Lemma (LS is satisfied, then (6.5) 
and (6.4 ) hold for n if they hold for n — 1. 



Proof. Consider ( 7.11 ) in order to prove that (3.2) holds for n. We have 



V fe (s - T ) = ^2 a kJ (s - T y 
i=i 

where the a^j are polynomials in e,, e$ and i = 0, ...^m^ (observe that we can 
consider the coefficients to be polynomials in ^„ ld , eo, e m<J and the if we wish). 
We thus get 



/ Vfc(g - r ) exp[-27(s - T )]ds = ^ / a fcj - 

•/ Tn „• 1 J TO 



(s - tq) 1 * exp[-27(s - T )]ds < 



< 



E 



23+iy+i 



We wish this expression to be less than or equal to 1, and since a k .j is polynomial in 
e$ , €i and ^ with each term containing at least one factor e, , we conclude that ( |6.3[ ) 
holds assuming that an inequality of the form ( [7.12 ) holds. By (2.6) and Lemma 
3.3, we have 



|P„,r(T,0)-P n;T (T O ,0)| < [ T [e- 2s \AP n \ + e 2P ^(Ql T + e- 2s \VQ n \ 2 )]ds < 
< A 7 /2 + 2 f Q\ exp[-2 7 (s - T )]ds, 

J Tn 



by Lemma 6.3. The remaining statement of the lemma follows. 



8. Convergence 



Consider the difference between (2.7) for n + 1 and n. We have 



Or. 



-It 



AQn+i + 2(P„ T Q n +i T — e 2t VP„ • VQ n +i 



2[Pn ,tQ n : t 



r VP n • VQ„) 
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(8.1) 



Qn+l(jo,-) = 
Qti+1,t(to,-) = 0. 



Since (6.4) is fulfilled for all n > and fl6.1| ) holds, Lemma |5.2| yields 
dE n k 

|c|=fe- 



(8.2) 

where 



< -2jE n . k + V / f n , a D a d T Q n dO, 



U,a = d T (e 2P "^D a d T Q n ) - V • (e 



2P„_i-2ty7 n a 



V£> Q Q„). 



We have 
(8-3) 

/n+i,o = e 2Pn [Qn+i,TT — e 2t AQ„ + i + 2(P„^ T Q„ + i :T — e 2t VP„ • VQ„+i)] = 

26 ^P n T Q n T 6 VP^ • V Qri) • 

The recursion formula is the same as before: if we have a — a + e;, then 

fn+l,a = dlf n +i^a ~ 29iP n / n + i^ — 

-2did r P n e 2P "D & d r Q n +i + 2e 2P "- 2T Vd i P„ • VD & Q n+1 . 

Lemma 8.1. Assume that conditions as in Lemma 7.4 are fulfilled, so that ft6.§ )- 
l \6..'\ ) are fulfilled for all n > 0. Let < k < nrid, then 

(8-4) E]ll xk < £ l n (lJ> k {T - r ) exp[-( 7 - 2A 7 )(r - r )], 

where Vk is a Q-dominated polynomial and £~n(m&x is given by (\6.{\). Furthermore 



.5) £ {\\e Pn D a d T Q n+1 \\ c{TdJ 

\a\—k 



V-D a Q n +l||c(T' i ,R)} < 



< £n(L*Qk+i(T - r ) exp[-( 7 - 2A 7 )(r - r )] 
on I , for k — 0, m^/2, where Qk is a Q-dominated polynomial. 



The proof of this statement is similar to the proof of Lemma 7.1. The proof of the 
following lemma is analogous to the proof of Lemma 7.2 

Lemma 8.2. // 

(8.6) ||e~ p "/, i+ i : a|| L 2 (T d iR) < £i/ 2 lax 7r fc exp[-( 7 - 2A 7 )(t - r )] 

for all a such that \a\ = k, where Hk is a Q-dominated polynomial independent of 



n, then an inequality of the form ( \8.4\ ) holds. 

Proof of Lemma \8.\ . Observe that for |a| > 1, 
(8.7) 

/«+!,«= A aua2 e 2P -{D^d T P n D a -d T Q n -e- 2T VD^P n -^D a -Q n ) + 

+ B Pu p 2 e 2P -{D^d T P n D^d T Q n+1 - e- 2T VD^P n ■ VD^Q n+1 ). 

/3i+/3 2 =a,|/3 1 |>l 
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In order to be able to apply Lemma S.2, we need to estimate the left hand side of 
(p!q). By the above observation (8.7), it is enough to estimate the expressions 



>D^d T P n D a *d r Q n \\ L * {T * 



VD^P n -VD^Q n \\ L2{Tit 



and 
(8.9) 



\e p "D^d T P n D^d T Q n+1 \\ L2{Td . m + \\e p »- 2T X7D^P n ■ VD^Q n+1 \\ L 2 



(T d ,R) ■ 



Concerning expressions of the form ( p.8[ ), we need only apply Lemma 5.1, Lemma 
7_+, (6J5) and the fact that one of |ai| and \a 2 \ has to be less than or equal to md/2 
in order to bound (B.5) by the right hand side of (8.6). Let us illustrate on the first 
term in (8.8) under the assumption that |ai| < md/2. We have 



\\e p "D^d T P n D^d T Q n \\ L 2 {T ^ 



< 



< 



\C(T d ,] 



\D a ^d T p n \ 



C(T d r 



D a *d T Q n \ 



L 2 (T d r 



< 



<Cexp[2A 7 (r-T )]f: 

If a = 0, then only terms of the form 
follows. 



1/2 v 

n,max' \oli 



exp[-7(r- r )]. 



are of relevance, so in that case, (3.4) 



Intermediate energies. Let us now prove (8.4) for < \a\ < m d /2 by induction 
on \a\. All we need to do is to prove that (8.9) is bounded by an expression as in 
the right hand side of (BJ3). Assume ( |S.4|) holds up to and includin g k < m d /2 — 1. 
Since k + 1 < md/2 and \[3i\ < \a\ = k + 1, we can use Lemma 5.1 in order to 
extract D^ 1 d T P n and e~ T VJ^f K in the sup norm. Since | /5s | < |a| — 1 < k, we 
can estimate what remains using (8.4) and the inductive hypothesis. By applying 
Lemma 8.2 we have thus proven 



4) for k < rrid/2. 



High energies. The argument now proceeds as in the proof of Lemma 7.1. By 



Lemma 6.1, D a P n is bounded in the sup norm by a polynomial if |a| < m,j/2, so 
that the result concerning the intermediate energies yields the conclusion that 

|l^ Q (e P ™<2n+l,r)||L 2 (T d ,K) 



can be bounded by the right hand side of (8.5) if |a| < Wd/2. A similar statement 



holds for e " T VQ n +i, and (B.5) follows for k = by Sobolev embedding. Assume 



now that (8^) holds up to and including k = m.d/2 + l < rn^— 1 and that ( |8.5| ) h olds 
for k up to and including I. For / = 0, we know this to be true. Consider now ( |S.9|) 
with /3i + (3 2 = a, \fix\ > 1 and \a\ = k + 1 = m d /2 + 1 + 1. If IftJ < m d /2, we can 
take out D^d T P n and e" T VD ft P„ in the sup norm by Lemma |5.l| and then apply 
the induction hypothesis to what remains since |/?2| < |a| — 1 < rrid/2 + l. If |/?2| < I, 
we can take out e p ^D^d T Q n+1 and e p "" T VO A Q n+ i in the sup norm, using the 



inductive hypothesis. What remains is bounded by £ n and (8^) follows for 
md/2 + 1 + 1 since one of the inequalities |/3i| < m d /2 and \j32\ < I must hold. In 
order to prove (|^) for I + 1, we proceed by Sobolev embedding as before. □ 

Finally, let us consider £ n ^. 
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Lemma 8.3. Assume that conditions as in Lemma 7.4 are fulfilled, so that ( ffi.^ - 
(\6.^ ) are fulfilled for all n > 0. We have 

< £ l JX* f V k (s - T ) exp[-2( 7 - 2A 7 )(s - r )]ds, 

J Tn 



(8.10) 



on I , for k = 0, m,d, where Vk is a Q-dominated polynomial. 



We need some preliminaries. Consider the difference of ( p.6[ ) for n + 1 and n. We 
have 

P„+i,rr - e- 2T AP n+1 = e 2P ™ (Qn+i, r - ^ T I VQ„+i | 2 )- 
-e 2P "- 1 (Qn,r " e~ 2r |VQ„| 2 ) = e 2P "(Q n+ i, r Q„+i,r - e- 2r VQ„+i • VQ„+i)+ 
+e 2P "(Q„, T Q n+ i, T - e- 2r VQ n • VQ„+i) + (e 2P " - e 2P ^)(Q 2 lT - e~ 2T |VQ„| 2 ) = 

= 9i + 92 + 93- 

Estimate 

(8.11) dS^k^y, r D *( gi+g2+g3)D * dr p n+ld0 

dT \ a \=k JTd 

Lemma 8.4. If 

(8.12) \\D a ( 9l + g 2 + «7 3 )IU»(r*,R) < ^max^kl exp[-2( 7 - 2A 7 )(r - r )] 



/or |a| < raj, where TT\ a i is a Q-dominated polynomial, then (8. It) follows. 
Proof. By Q , 

< 



a I— A' 



which can be integrated to ( 8.1C ) given the assumptions of the lemma. □ 



Proof of Lemma 8.S. Consider the contribution of the first term in g\ to D a g\. If 



more than md/2 derivatives hit one of P n , Q n +i.r or Qn+i,r, then everything else 



can be taken out in the sup norm, and we obtain an estimate of the form (S.12) 
In fact the estimate is a bit better, but we will need the form ( |B.12 ) for the other 



terms. The argument for the contribution of the second term in gi is similar. The 
expression D a g2 can be controlled by similar arguments, but we loose one factor 
exp[— 2Aj(t — To)] in decay due to the fact that we have to compensate that we in 
some situations have e p ™ were we would prefer to have e p ™ _1 . 

Consider 

D « [(e 2P„_ e 2P„_ 1)Q y. 

This expression is a linear combination of the following types of terms 

(e 2P "^ Ql P„ • --D^P n - e ap »- 1 D Q1 P w _i ■ ■ ■ D a ' P n ^)D a '+^d T Q n D a '+-d T Q n 

where ai+...+a;+a = ol. These terms can in turn be written as a linear combination 
of terms such as 

(8.13) e 2P - 1 D a 'P n ■ ■ ■ D^P n ■ ■ ■ D a <P n ^D a ^d T Q n D a ^d T Q n 
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and 
(8.14) 



(e 2P " - e 



2P„_i \ noi 



)D ai P n ■ ■ ■ D ai P n D a ^d T Q n D ai + 2 d T Q r , 



Terms of the form ( p. 13 ) can be handled using Lemma 6.1 and 7.1 , since at most 
one | ay | can be bigger than rrid/2. For terms of type ( |8.14 ), we estimate 



| (g2 p„ _ e 2Pn- 1)D a iPn . . . D a >P n D a ^d T Q n D a ^-d T Q n \ < 
< 2|P„ - P n _i| exp[4A 7 (r - r )]e 2P "- 1 |D Q1 P„ • ■ • D°" P n D a '^d T Q n D a ^d T Q n \ 



which we can deal with, using (6.11). The contribution from the second term in 173 
can be estimated similarly. □ 



9. Conclusions 

Theorem 9.1. Let < 7 < 1/4, t > and define m d = 2[d/2] + 2. There are 
constants c%,d) i = 1, 2, and integers li t d, i — 1,2 depending on the dimension such 
that the following holds. If d — 1, let £2 and To be such that 

(9.1) (l + e 2 )e- T0 <c M7 . 

If d > 1, let v md / 2 +i, e mj /2+2 and r be such that 

(9-2) + tm d /2+2 + l]e~ 2r ° < Cl )d 7. 

Specify the remaining €k and Vk, k = 0, freely. Assume furthermore that 

(9.3) e fc < c 2)d (l + e + e md + ^J-^V 2 '" 

/or A: = 0, ...,md- Then every quadruple of functions (poiPii <7o, 9i) satisfying 

(9.4) 27 < pi < 1 - 27 



and {6.1), yield upon solving \2.\ ) and (2J<) smooth solutions on [to, 00) with the 
following properties. For all non-negative integers fc, there are polynomials 'ELi.k, 
i = 1, 6 in t — t , and v, w,q,r € C°°(T d , K), where 

(9.5) 

on T d such that 



0<7<f<l~7<l 



(9.6) 
(9.7) 

(9.8) 
(9.9) 

(9.10) 

and 

(9.11) 



IIP- - f||cHT d ,R) < s i,fcexp[-27(r-ro)], 
ll-P - Hlc fe (T<*,K) < s 2,fe exp[-27(T - t )], 
J] \\e 2p D a d T Q\\c(T*X) < S 3 , fc , 

|a|<fc 

2 \\e 2p [D a Q-D a q]\\ c{Tdm < ~ 4 , fc , 

a|<fe 

e 2p Qr - r|| C fc (T d :K ) < 3 5)fe exp[-27(T - r ) 



|e 2p (Q-g) + — 



< S 6jfe exp[-27(T - t )] 
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for all t £ [ tq,oo) ; where p = v ■ (r — To) + w. Finally, assume v G C°°(T d ,R) 
satisfies (9.1), and that 



(9.12) \\e P D a d T Q\\ 

i«i<"i<j+i 



C m d + 1 (T d , 



< Ce" 



w/iere e > and P, Q G C°°(R x T d ,R). Let 

ej t (r)=^ /a (P,Q,r) ) e' fe (r) = ^ /a (P,r), i/ fc (r) = ^ k /2 (P, r) 



and 7' = 7/4. Then for Tq big enough, (9.1)-(9.4) will be satisfied with e k replaced 
by b'^Tq), tq replaced with Tq, 7 replaced with 7' etc. 



Remark. The last part of the theorem is intended to illustrate the fact that there is 
in some sense an equivalence between the conditions ( 9. 1 )-( 9.4 ) and the asymptotics 
(PI)-(lPll). 

Let us briefly comment on the conditions before we turn to the proof. The idea of 
the argument is to see to it that P T is always bounded away from below and 1 
above. As can be seen by Lemma |6.3| , the conditions fl9.l| ) and (9.2) are there to 
ensure that the term e _2r AP in (2.1) does not push P T out of this interval. The 
condition fl9.3Q is then there to ensure that the remaining terms in (2.1) do not 
push us away from the desired region. 



Proof. Observe first that (6.3) and ( p.4| ) hold for n = due to Lemma 6.4 and the 
assumption that ( |9.l[ ) or (9.2) hold. By Lem ma \IA\ an estimate of the form (9.3), 
together with estimates of the form ( |9.l| ) or (|9.2|) imply that (6.3) and (6.4) hold 
for all n > 0. By (|8.1C), an argument similar to the proof of Lemma 7.4 yields that 



a condition of the type (|9.3|) implies 



(9.13) 



£1/2 < 1 £1/2 

n+l,max — 2 n ' max " 



Let T G / and let us consider the convergence on [to,T] x T d . By (3.11), we 
conclude that P n is a Cauchy sequence in sup norm. Observe that under these 
circumstances, factors of the type e 2Pn and e _2r are of no importance, since we are 
considering a finite time interval, and since the sequence P„ is uniformly bounded 
on this finite time interval. By ( |8.5D , ( |9.13 ) and the equations, we conclude that 



D a 8 l T P n and D a d l T Q n 

are Cauchy sequences in C([tq, T] x T d , K) for 1 < |a| + 1 < md/2 + 1. The conver- 
gence of Q n follows from the convergence of Q n ,r, the finiteness of the time interval 
and the fact that the iterates coincide at to. In particular, t he it eration yields 
C 2 solutions to the equations for t G [to, 00). By Proposition 3J, the solutions 
will be smooth if the initial values consist of smooth functions (the transformation 
t = —e~ T yields an equation of the right form). Furthermore, if we let 



and 



E k {r) =E k (P,Q,r) 



£ k (r) = £k(P,r), 



then these expressions will satisfy the estimates (7.1), (J7. 2|) and (6.3) 
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Lemma 9.1. There are ■polynomials Vk and Qk and constants Ck such that 
eI /2 (t) < V k (T - r ) cxp[- 7 (T - r )], 

{\\e P D a d T Q\\ c{TdM) + \\e p - T VD a Q\\ C{Tim } < 

\a\—k 

< Qk+i{r - t ) exp[- 7 (r - t )] 



(9.14) 
(9.15) 



and 
(9.16) 



£k < Ck < oo 



on I, for all non-negative integers k. 



Proof. By our construction we know the statement concerni ng (9.14 ) and ( 9.16 ) 
to be true for k = 0,...,rrid and the statement concerning (9.15) to be true for 
k = 0, ...,mrf/2. We want to prove the statement of the lemma by an inductive 
argument. Assume it to be true up to and including k > for (9.14) and (9.16) 
and to k — rrid/2 for (9.15). Let us introduce 



\WD a P\ 2 d9 



\ol\— k ' 



and 



Qk = exp[-7(r - T )]F k - 
Our primary goal is to prove the following two inequalities 



(9.17) 

and 

(9.18) 



dE, 



fe+i 



< -2 7 £ fc+1 + [y k+ i,i + ^+i,2^{i] exp[- 7 (r - 7to)]Ej£i 



,,1/2 



dr 
dSk+i 



< [Zk+is + Zk+iaE]/^ + Zk+x^g 1 ^] exp[- 7 (r - tq)]£^ x 



dr 



l/2n 



,1/2 



where 34+1, i> 34+1,2, Z k +is, Zk+1,2 and Z k +i^ are polynomials in t-t . Observe 
that we also have 



(9.19) 
We have 
(9.20) 

where 



dS, 



< -jGk + C k exph 7 (r - r )/2]^ /2 £ fc 1 / 2 1 . 



dr 

dEk+i 
dr 



< -2jE k+1 



E 

|a|=fc+l 



f a D a d T Qd6 7 



'I 'd 



/« = UP, Q) = d T {e 2F D a d T Q) - V • ( 



„2P-2t 



VD"Q) 



and if a — a + e; , 
(9.21) 



f a =d l f & ~28 l Pf & - 
-2did T Pe 2P D & d T Q + 2VdiP ■ e 2p - 2r WD & Q. 



Let us first prove (|9.17| ). Considering (|9.17| ) and (|9.20| ), it is enough to prove the 
estimate 



(9.22) 



'fa\\i?(T* ,K) < [IIfe + i,i + Tlk+i,2 £ l+i] exp[- 7 (r - t q )} 
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when \ot\ = k + 1, where Ilfc+^i and Tlk+1,2 are polynomials. Since /q is zero, we 
inductively get the conclusion that f a will contain two types of terms: 



(9.23) 

and 

(9.24) 

where |«x| > 1 an d ot\ 
we wish to estimate 



On 



„2P 



{D ai d T P)(D a2 d T Q) 



-ol-2 — a. Consider terms of type ( 9.23| ). Considering (9.22) 



\e p {D a -d T P){D^d T Q)\\ L , (Td ^. 



If la x] < k — 1714/2, then we can estimate D ai d T P in the sup norm using Lemma 
5.1 and the inductive hypothesis conc erning ( |9.16| ). Since \a^\ < \a\ — 1 < k, the 
inductive hypothesis concerning ( |9.14 ) yields the conclusion that 

\\e p (D a 'd T P)(D^d T Q)\\ L 2 (Td , R) < C k+1 V {a2l eM-l{r-r )]. 

If |cK2 1 < k — rrid/2, then we can take out e p D a2 d T Q in the sup norm, using the 
inductive assumption concerning (9.15), in order to obtain 

\\e p (D^d T P)(D^d T Q)\\ L 2 (TdM) < ^Q| Q2 | +1 exp[- 7 (T-r )]£ | 1 Q / i 2 | . 

Regardless of whether |ai| = k+ 1 or not, we get an estimate that fits into ( |9.22 ). 
As we cannot have |aa| > k — rrid /2 an d |ai| > k — m^/2 at the same time, we 
have dealt with terms of the form (9.23). Terms of the form (3.24) can be handled 
similarly. Equation (9.17) follows. 

As far as ( |9.18 ) is concerned, we have 
(9.25) d£k+1 



dr 



< 



\a\=k+l 



D a [e 2P (Q 2 T - e- 2T \VQ\ 2 )}D a d T Pd0. 



Considering (9.18), it is thus enough to prove 



(9.26) 

for \a\ = k 
(9.27) 



a\„2P 



-2t 



|vqi 5 



< [n fc+ i,i + nfc+1,2^ + n fe+ i, 3 ^ /z ] exp[- 7 (T - r )] 
1 and some polynomials Ilfc+i 1, IIfc +1 2 and Hk+i,3- Consider 



|i 2 (T<i,: 

1/2, 



< 



D a (e 



IP, 



k+1 

= Y, J2 C au ..., ai+2 e 2P D ai P---D "PD°"+ 1 d T QD a '+ 2 d T Q. 

1=0 a 1 + ...+a i+2 =a 

Observe that we can bound D°P by a polynomial in supremum norm if \/3\ < 
k — rrid/2 and in L 2 norm if |/3| < k using the induction hypothesis concerning 
(3.16) and Lemma 5.1. We can also control e p D l3 d T Q in the sup norm, using 
(3.15) and the inductive assumption, if \/3\ < k — and in the L 2 norm if 



\/3\ < k. Consider a term in the sum (9.27). At most one |a,| can be greater than 
k — nid/2. If all on satisfy \ol\\ < k, then we get a bound 

\\e 2P D ai P ■ ■ ■ D a 'PD a ^d T QD a ^d T Q\\ L 2 {Td , R) < U k+1 exp[-2 7 (r - r )] 

where Ilfc+i is a polynomial. If = k + 1, and i < I, i.e. if k + 1 derivatives hit 
one P, then 

\\e 2P D ai P- ■ ■ D a ' PD a ^d T QD a '+ 2 d T Q\\ L 2 (Td , m < 
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< n fe+1 exp[-2 7 (r - T )]\\D a P\\ LHTd , m < U' k+1 exp[-3 7 (r - r )/2]^/ 2 . 
Finally, if one \eti\ = k + 1, and i > I, i.e. if all the derivatives hit one d T Q, then 

\\e 2P D ai P- ■ ■ D ai PD a ^d T QD a '+ 2 d T Q\\ L 2 (Tdm < 

< n fc+1 exp[- 7 (r - r )]||e p £) a 9 T Q|| L 2 j:TI i im) < W k+1 exp[- 7 (r - t q )]E\^ v 

As the argument for D a (e 2P ~ 2T | VP| 2 ) is similar ( |9.26| ), and thereby fl9.18| ), follows. 
Let 

riu = £k+i + Ek+i + Gk- 
The estimates (|9.17| )-( |9.19D imply 

(9.28) ^± < W kA exp[- 7 (r - r )/2] r^/ 2 + W fe , 2 exp[- 7 (r - r )/2]W fc , 

CtT 

where Wk,i and Wk,2 are polynomials. This inequality implies that TLk is bounded 
on /. In fact, since a 1/2 < |(1 + a) < 1 + a for a > 0, ( |9.2g| ) yields 

d(l + Wfc) 



< {W k ,i exp[- 7 (r - t ) /2] + W k , 2 exp[- 7 (r - r )/2]}(l + H k ) 



whence 1 + TLk is bounded. Thus £k+i is bounded, yielding (9.16) for k + 1, which 
when inserted in (9.17) implies (9.14) for fc + 1, which, together with the induction 
hypothesis yields ( |9.15| ) for k + 1 — md/2 by Sobolev embedding. □ 

Observe that using ( 9.16| ) and Lemma 5.1, we can conclude that 

\D a P\ < K a 

on / x T d for all a, whe re lZ a are first degree polynomials in r — tq. Using this 
together with (2J) and ( 9.15 ), we conclude that 

\D a d 2 T P\ < T Q exp[-2 7 (r - r )] 

on / x T d for all a, where F Q are polynomials in r — tq. This implies that d T P 
converges in any C k (T d ,M.) norm as r — > oo to a smooth function u on T d , and 
that ( |9.6| ) holds. Equation ( |9.6[ ) in its turn implies that 

P(r, •) - v ■ (t - t ) w 

in any C k (T d ,M.), and that we have (|9~7|). Observe now that by (2.2) we have 

d T (e 2P d T Q) = V • (e 2p - 2r VQ). 

Thus 

|a rJ D Q (e 2P Q T )| = |V • D a (e 2p - 2T VQ)\ < ir a exp[-2 7 (r - r )] 



combining (9.15) with the fact that P—r < C — 7 (t — tq). Here and below, n a will 
denote a polynomial in t ~ tq. This inequality can be integrated in order to yield 
the conclusion that there is an r £ C°°(T d ,M.) such that 

(9.29) \D a {e 2P Q T -r)\ < tt q exp[-2 7 (r - r )]. 

Let now 

P = v ■ (t - r ) + to, 

and observe that we have (3.7). We would now like to replace P in ( 9.29| ) with p. 
Consider for that reason 

D a (e 2p Q T -r)= D a [(e 2 P - e 2P )Q T ] + D a (e 2P Q T - r). 
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The second term we can bound using ( |9.29 ). Let us consider the first term, 
D a [{e 2p -e 2P )Q T ] = £ C ai . a2 D a - {e 2p - 2P - l)D a *(e 2P Q T ). 

The second factor on the right hand side is bounded by a constant for r > To by 
the above, and so the factor of interest is 

D ai (e 2p - 2P - 1). 

Consider first the case a± =0. If r is great enough, then \2p — 2P\ < 1, so that 
|e 2 P -2P _ i| = | g (2P ~^ P) " | < 2\p - P\ exp[|2p - 2P\] < 2e\p - P\ < 



< 7T exp[-27(r - r )], 

where we have used fl9.7|) in the last inequality. If |ai| > 1 we get similar estimates 
for less complicated reasons. To sum up, 

\D a {e 2p Q T ~r)\ < ix a exp[-2 7 (r - To)], 

proving fl9.10|) . We conclude that holds, and therefore there is a q e C°°(T d , E) 
such that (3.£) holds. Let us compute 



e 2p (Q - q) = e a '( 



2p/ 



Q T {s)ds) = - 

e 2p ^- 2p ^(e 2 P^Q T (s) -r)ds 
Observe that r is independent of s and compute 

/>oo 

3 2p(r)-2p( S ) ds = / e -2v(s-T) ds 



2p(r) 



} T (s)ds 



,Mr)-2 P (s) rds _ 



2v 



We thus conclude that 



D a {e 2p (Q -q) + ^-) = -D a Qf°° e^-^V^Q^s) - r)ds^ 



/ />OC 

X! C ai , a2 ( D ai [e 2p ^- 2p ^}D a2 {e 2p ^Q T 



(s) — r)ds 



ai-\-a2—a 

In order to estimate this expression, consider 



| / D ai [e 2p ^- 2p ^]D a2 (e 2p ^Q T (s)-r)ds\ < 



< 



[D ai e 2p{T) - 2p(s) ] 2 ds 



1/2 



[D a2 (e 2p< - s) Q T (s)-r)] 2 ds 



1/2 



The integrand in the first factor can be bounded by a polynomial in s — t multiplied 
by exp[— Av ■ (s — r)]. In consequence the first factor is bounded by a constant. The 
integrand appearing in the second factor can be bounded by 

tt Q2 exp[-47(s - r )] 
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where polynomial in s — to- Adding up these observations, we get the 

conclusion 

\D a (e 2 ?(Q -q) + ^-)| < ir a exp[-2 7 (r - r )], 



proving (3.11) 



Let us now prove that solutions with the desired asymptotics satisfy the initial 



conditions at late enough times. Due to (9.12) there is a w £ C" ld+1 (T d ,M.) such 
that 

(9.30) \\P(t, ■) - v ■ (r - r ) - HIc-^+H^R) - « 

as t — > oo. As a consequence, the e' fc (r) are bounded for the e ntire futur e and the 
v' k {r) do not grow faster than linearly. Furthermore, by (|9~30| ) and (fyj), 



exp(— w • r — er) 

a|<m d + l 

Thus there is a 9 G C md+1 (T d ,R) such that 

^ ||e p ^ Q (Q- g )|| c(T£iiR) <Ce- er . 

\a\<m d + l 

Since P < C + (1 — 7)t and 7 > 0, we get the conclusion that e' fc (r) decays to zero 
exponentially in time. The last statement of the theorem follows. □ 



10. Curvature blow up 



Let us make some observations concerning the geometry of the metric (1.1) given 
the conclusions of the previous section. 

Proposition 10.1. Consider a metric of the form fli.^ ). Assuming P and Q have 



the asymptotic behaviour obtained as a conclusion in Theorem 9.1, we have 



lim inf \(R a 3jsR°' 0lS )(T,0)\ =00. 
Proof. We proceed as in [|| . Consider the orthonormal basis given by 

eo = e A/4+3r/4^ ) ei = e A/4-r/4 %) ^ = ^-PI^ e 3 = e T/2+P/2 {d S - Qd a ). 

There is a natural scaling given by 

4> = e -V4-3r/4_ 

Only few terms in (fpRa^s, where we assume that the indexes are with respect 
to the orthonormal basis mentioned above, are non-negligible. This makes the 
computation of the Kretschmann scalar manageable and in fact, 

converges to a smooth non-zero function on S . Since A equals v 2 (6)t up to some- 
thing bounded, we get the conclusion of the proposition. □ 

There is also an elementary proof of causal geodesic incompleteness in our situation. 
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Proposition 10.2. Consider an inextendible causal geodesic 7 : (s_,s_|_) — > M, 
where M = M x T 3 wii/i a metric of the form (1.1) where P and Q have the 
asymptotic behaviour obtained as a conclusion in Theorem 9.1. Assume furthermore 
that < 7 '(s), 9r|-y(s) >< 0- 27^71 7 is future incomplete and 

lim |(i? Q ^i? Q ^)[ 7 ( S )]| =00. 

S— »S-|- — 

Proof. Let the basis e M be as in the proof of the previous proposition and define 

fo = - < 7', e o| 7 >■. /fc 7 ',e fc | 7 > 

for k — 1,2,3. Observe that Ylfk — fo ^ ue ^° causality. Furthermore, if we let 
the t component of 7 be denoted by 70, then djo/ds > 0. Thus, if 70 is bounded 
from above, it converges to a finite value as s — ► s+— . By the causality of the curve 
and compactness of the spatial slices, we then get the conclusion that 7 converges. 
Thus 7 is continuously extendible, leading to the conclusion that it is extendible 
considered as a geodesic. This contradicts our assumptions, and we conclude that 
7 o(s) -^00 as s — > s+ — . Consider 

df 

= ~ < l', V 7 'e >= - ^ ffifv < e M ,V e „e > 07. 

Let <p be as in the previous proposition. Using the assumptions concerning the 
asymptotics one sees that 



and 



< ei, V ei e >= -~(A T - 1), <f) < e 2 , V e2 e >= ~(1 - P T ) 



< e 3 ,V e3 e >= ~(1 + Pr) 



and that all other elements of the matrix <f> < e„, V ei/ eo > converge to zero expo- 
nentially with t. Letting 

Ok = 4> 7 < efe, V efc e > 07, 

we have 

*=-^tE flh + 4' o 7( 5i/ 2 > 

k 

>4>° 7/ot( a t °7 - 1) + ip 7^i/oi 

if s is close enough to s + , where ip = l/(f> and S%(s) — ► as s — ► s + — (observe that 
6*i > and 2 , #3 < if s is close enough to s + ). Compute 
dip o 7 9-0 d 7o 9-0 d 7 i 
ds 9r ^ ds 89 ^ ds 
where 71 is the ^-coordinate of 7 (observe that even though this is not well defined, 
the derivative is). However, 

djo , , dli 



so that 



^°lfo, -7 - =exp[A 07/4-70/4]/! 
ds ds 



ON GOWDY VACUUM SPACETIMES 



27 



where ^(s) — > as s — > s + — . Letting g = fo ■ ip o 7, we get 

-j- = -7-^ 7 + /o j- > .g [j(A r o 7 - 1) + Sx + — (A T o 7 + 3) + 6 2 }. 

Thus there is an si such that for s > si 

d 9 ( \ ^ 1 2 / \ 
^ (S) ^2 3(S) ' 

since A T is bounded from below by a positive constant for large r. We thus get the 
conclusion that the geodesic is future incomplete. Since 70 (s) — > 00 as s — ► s+— , 
the statement concerning curvature blow up follows from the previous proposition. 
□ 



11. Conclusions and observations 
This paper provides a proof of the statement that an open set of initial data yields 



asymptotics of the form (1.6) and (1.7). The fact that we have a condition on initial 
data will hopefully be useful when trying to make further progress in understanding 
Gowdy spacetimes. However, the method of proof relies heavily on the fact that the 
non-linear terms in the equations die out exponentially. This will not be the case in 
general and different methods will have to be developed. Some intuition for what 
should happen has been developed in but how to make these ideas rigorous is 
far from clear. Since approaching the general problem is difficult, it is natural to try 
to find an easier problem which can be used as a model for some of the dynamics. 
One way of obtaining a model problem is to carry out intuitive arguments similar 
to those in 0]. For example, if P T > 1, e p Q T should be negligi ble a nd Qg should 



be a time independent function. Inserting these assumptions in (1.2) one obtains 

r TT — e i^gg — — e y e . 
Replacing Qg with 1 and calling P — r P, one gets the equation 
(11.1) Prr - e- 2T Pgg = ~e 2P . 

This equation would, in the terminology of H, model a bounce. Considering an 



arbitrary solution to (11.1), P T should converge exponentially to a smooth negative 
function on S 1 . One question of interest would then be to find out how fast this 
process of convergence occurs. This is an example of a model problem which is 
certainly easier to handle, but which still contains some of the dynamics of the 
real situation. Since most of the results concerning the asymptotic behaviour of 
Einstein's equations in the spatially inhomogeneous case make assumptions that 
exclude the possibility of the non-linear terms being of any importance, the study 
of even a simple model problem where this is not the case should be of interest. 
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